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A piecewise continuous map for modeling bursting and spiking behaviour of isolated neuron is 
proposed. The map was created from phenomenological viewpoint. The map demonstrates oscil- 
lations, which are qualitatively similar to oscillations generating by Rose-Hindmarsh model. The 
synchronization in small ensembles of the maps is investigated. It is considered the different number 
of elements in the ensemble and different connectivity topologies. 
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I. INTRODUCTION 

Investigations of neuron ensembles (Central Pattern 
Generators, CPG) from the nonlinear dynamics view- 
point attract attention of many physicists last years. It is 
known some successful experiments with biological neu- 
rons (see, for example, [1]). Now exist a lot of math- 
ematical models describing neuron behavior on the ba- 
sis of ordinary differential equations [2]. Even simple 
models such as well-known FitzHugh-Nagumo and Rose- 
Hindmarsh systems, which may describe some important 
facts of neuron dynamics, need a lot of resources in nu- 
merical experiments for rather large ensembles. In this 
case the models with discrete time (maps) are much more 
appropriate. 

It is known a set of such kind maps (i) the model with 
variable taking the discrete set of values (finite automata) 
[3,4], (ii) the map taken on the flow of differential equa- 
tions system [5] , (iii) learning globally coupled excitable 
map system [6], (iv) recently proposed two-dimensional 
map [7]. However, the variety of models and methods for 
the neuron dynamics descriptions testifies that nonex- 
ists the common approach for creation of maps defining 
the isolated neuron behavior as well as for creation of 
the model neuron ensembles with pregiven topology and 
coupling kind. 

In this paper is proposed a map with variable that 
changes continuously in specified range. This map 
demonstrates behavior that is qualitatively similar to 
the real neuron dynamics. It will be shown below that 
proposed map describes the complete synchronization in 
neuron ensembles. It will be also observed dependence 
the synchronization degree of the network configuration. 
Another results that will be discussed connected with the 
influence of external force on studied ensembles. 

II. THE SINGLE NEURON MODEL 

The map was constructed on the assumption of phe- 
nomenological conceptions. The basic idea is founded on 



the fact that in neuron behavior at time series the three 
regions (rest, burst and spike regions) could be identified. 
This relative division is schematically shown in Fig. 1. 

The map (x — > x') consists of two piecewise continuous 
functions connected by transition conditions. The system 
dynamics is described with four state variables: x G [0, 1], 
d = { — 1,1}, si,s 2 = {0,1}- The main "observable" 
variable x qualitatively corresponds to membrane poten- 
tial of the neuron, the variable d is used to select one 
of two branches of function and the variables s\,S2 are 
"switches" , which define the conditions for burst ending. 
In mathematical notation the map are written in the fol- 
lowing way. 
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FIG. 1. The fragment of time series of isolated biological 
neuron ("1" is the rest state region, "2" is the burst and "3" 
is the spikes). 

If d = 1 (the increasing of the x value) 

!ai arctan(fcia;), if x G [0, a — Si) , 
2a — x, if x G [a — Si, a) , (1) 

71 (x — a) + a, if x G [a, c) , 

The transition condition: d = — 1 if x G [c, 1] . 
If d = — 1 (the decreasing of the x value) 

!i (x — a) + a, if x € [a + 62, 1], 
2a — x, if x £ [a, a + S2) , S1S2 — I (2) 

arctan^ir), if x G [£3, a) . 
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Additional terms: d = 1 if x G [a, a + 5 2 ) and sis 2 = 0, 



si 



1, if a; G [c, fti], 
0, if a; G [0, a], 



S2 



l, if xe[h 2 ,i], 
0, if x G [0, a], 



d= 1 if x G [0,(5 3 ). 

Here a, fci, fe, 71, 72, <$2, <$3 are some con- 
stant parameters, and others coefficients are determined 
from continuity condition (c = (1 — a)/ji + a, a\ — 
a I arctan(fcia), a 2 — a/ arctan^a)). Parameter h 2 is 
also constant and close to unity, and hi should depend 
on the value of c, that is hi = c + Ah, where Ah is the 
relatively small parameter. 
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FIG. 2. The iteration diagram for the map, describing by 
(1), (2) transition conditions. 

For better understanding how this map works the iter- 
ation diagram is shown in the Fig. 2. Fig. 3 demonstrates 
the time series for the isolated model neuron described 
by the map for various values of control parameters. 



III. GENERALIZATION OF THE MODEL TO 
COUPLED NEURONS. TWO COUPLED 
NEURONS 

Investigation of single neuron behavior has a little 
practical sense. Much more interesting might be model 
of neuron ensemble. Therefore this piecewise continu- 
ous map was generalized to coupled neuron systems. In 
case of the ensemble with N elements with state vari- 
ables (x^jX^, ■■■,Xn) the influence of neighbours to j-th 
neuron at (n + 1) discrete time is given by adding the 



following term to x 3 n+1 variable 
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T 



N 

i = 1 



(3) 



where £y, = 1,...,N) is the connection weight be- 
tween i-th and j-th. elements, Lj is the number of neigh- 
bours of j-th neuron, 6 is the Heaviside step function 
(authors presume that exists a threshold of interaction). 



The first investigated small ensemble was constructed 
of two coupled neurons. The main interest of the system 
research was consisted in revealing complete synchroniza- 
tion. The degree of synchronization was determined as 
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(4) 



n=l 



where M is the number of iterations for a discrete time 
average. 
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FIG. 3. Characteristic time series for isolated model neuron 
for 82 = 0.001, hi = 0.88, 71 = 1.4, 72 = 1.75 (a); S 2 = 0.001, 
hi = 0.88, 71 = 1.402, 72 = 1.75 (b); 8 2 = 0.0001, hi = 0.92, 
71 = 1.3, 72 = 1.3 (c). The values of the other parameters are 
a = 0.3, fci = 0.9, k 2 = 1.0, h 2 = 0.95, 5i = 0.01, S 3 = 0.001. 

Value A = corresponds to the highest degree of syn- 
chronization, A > means that the synchronization is 
absent or synchronization is incomplete. In this research 
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the degree of synchronization (4) on the plane of param- 
eters (e, 71) is calculated. The choice of parameter e as a 
control parameter is natural because it's value determines 
the strength of connection. The value of parameter 71 (as 
well as 72) influences considerably to the system behav- 
ior, so it is chosen as second control parameter. The pa- 
rameter planes are plotted in grayscale, namely synchro- 
nization (for value A = 0) region marked with the white 
color, the positive value of A is subjected to the simple 
rule: the larger value, the darker point. Such parame- 
ter plane for two coupled neurons is presented in Fig. 4. 
The behavior of this system in some points, marked on 
the parameter plane with letters, is shown in Fig. 5. 
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FIG. 4. The parameter plane, plotted for two coupled 
model neurons. The letters correspond to time series, pre- 
sented in Fig. 5. The values of the other parameters are 
a = 0.3, fci = 0.9, k 2 = 1.0, 72 = 1.75, 61 = 0.01, 8 2 = 0.001, 
S 3 = 0.001, h 2 = 0.95. 

It was revealed that the system dynamics depends on 
the initial conditions therefore the question about typi- 
calness of presented parameter plane occurs. This prob- 
lem was studied in detail. The view of the attraction 
basins, plotted in different points of the parameter plane, 
and the parameter planes obtained for various initial con- 
ditions shows that the general features of the planes are 
invariable, so this question will be not discussed below. 

IV. SYSTEMS WITH SEVERAL ELEMENTS 

In the research is observed a large quantity of neuron 
ensembles with different numbers of elements and vari- 
ous spatial configurations. Some revealed regularities for 
ensemble of seven coupled neurons is presented here. 

In case of several elements the degree of synchroniza- 
tion can't be defined with (4), therefore in numerical ex- 
periments is used the following approximate relationship 
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FIG. 5. Characteristic time series plotted for two coupled 
neurons in the points, marked on the parameter plane in fig. 4 
with e = 0.50, 71 = 1.40 (complete synchronization) (a); 
e = 0.95, 71 = 1.70 (b); e = 0.24, 71 = 1.427 (c). 

Fig. 6 presents the parameter planes that are gotten 
for systems schematically shown at the insets. As it can 
be seen, the areas with synchronous behaviour become 
larger when additional connections are inserted. Note, 
that although this tendency is clearly visible from Fig. 6, 
the synchronous behaviour is observed in another areas of 
the parameter plane (see Fig. 6b and Fig. 6c). Comperi- 
son of Fig. 6b and Fig. 6c shows that synchronization 
degree is decreased in the region e s [0.8; 1.0] for all 
values 71. This result remains correct for systems with 
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another numbers of elements. 

It was shown that the synchronization degree de- 
pends on the fast motion (the spikes region) and prac- 
tically is independent of slow motions. This fact was 
revealed when the nonlinear branches (sec Fig. 2, region 
x n € [0; a}) of the map functions (1), (2) were approxi- 
mated with piecewise linear functions. There is the visi- 
ble changes of the time series, but the view of parameter 
planes remains without any considerable changes. 
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FIG. 6. The planes of parameter for ensembles with seven 
elements and different spatial configurations (shown at the 
insets) obtained from ring-type system by adding of connec- 
tions. 



V. THE INFLUENCE OF THE EXTERNAL 
ACTION 



shown that the proposed map qualitatively simulates real 
neuron behavior and describes some synchronization phe- 
nomena that is observed in neuron ensembles, connected 
via electrical synapse. The assumption that just fast mo- 
tion is responsible for the synchronization phenomena is 
prior and should be confirmed. The results related to a 
possibility of low-amplitude external field to increase the 
degree of synchronization in small model neuron ensem- 
ble may be useful from the practical viewpoint but only 
in case that they will get an experimental verification. 
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FIG. 7. The ring-type system with seven elements under 
the low-amplitude spatially uniform periodic (a) and chaotic 
(b) fields £ where a — 0.0003. The periodic field is 

generated by y„ — sin (2nn/T) with T = 500 and chaotic one 
by the logistic map y„+i = 1 — \y„ with parameter A = 1.99. 

The way on the creation of the piecewise map from phc- 
nomenological viewpoint (from view of the system time 
series), applied in the work, could be used for other sys- 
nems. 

This work was supported by RFBR under grants 02- 
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and Ministry of Education of Russian Federation under 
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The investigations of the synchronization also involved 
the possibility to influence on this process with an exter- 
nal signal. It was revealed that an external low-amplitude 
spatially uniform field (both periodical and chaotic) may 
increase the degree of synchronization in the ensemble. 
As the example of this fact the parameter planes are pre- 
sented for ring-type system with seven neurons under an 
external force (Fig. 7). 

As one can see from comparison of the parameter plane 
in Fig. 7 and Fig. 6a, the region of synchronous behaviour 
for the model neuron ensemble under the external force 
becames larger in relation to the region for unperturbated 
system. 

VI. CONCLUSION 

In this work the phenomenological map describing 
some aspects of neuron dynamics is presented. It is 
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